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Abstract — In this paper, we consider spatial-division multiple- 
access (SDMA) systems with one base station with multiple 
antennae and a number of single antenna mobiles under noisy 
limited CSIT feedback. We propose a robust noisy limited 
feedback design for SDMA systems. The solution consists of 
a real-time robust SDMA precoding, user selection and rate 
adaptation as well as an offline feedback index assignment 
algorithm. The index assignment problem is cast into a Traveling 
Sales Man problem (TSP). Based on the specific structure of the 
feedback constellation and the precoder, we derive a low complex 
but asymptotically optimal solution. Simulation results show that 
the proposed framework has significant goodput gain compared 
to the traditional naive designs under noisy limited feedback 
channel. Furthermore, we show that the average system goodput 
scales as £>( ";^'-_7' (C/b - log^ (jV„))) and ©(nr -loga P) in the 
interference limited regime (Cft < (nr — 1) logj P + logj A''„) 
and noise-limited regime respectively. Hence, despite the noisy 
feedback channel, the average SDMA system goodput grows with 
the number of feedback bits in the interference limited regime 
while in noise limited regime increases linearly with the number 
of transmit antenna and the forward channel SNR (logj P). 

I. Introduction 

It is widely known that spatial-division multiple- 
access(SDMA) is an important technique to enhance 
the throughput of multi-user wireless systems due to spatial 
multiplexing. However, SDMA system requires channel state 
information at transmitter(CSIT). In FDD systems, only a 
limited number of bits(e.g. 6 bits for WiMAX [3|) can be 
allocated to carry the CSIT feedback, namely the limited 
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feedback. In the authors consider multiuser MISO system 
with a limited total feedback bits constraint and proposed 
a codebook design algorithm and a CSIT decomposition 
algorithm. The authors of ^ analyzed the asymptotic 
performance of a per user unitary and rate control design. 
Other works like ||4l-||5l studied transmit beamforming using 
different criteria and methods. However, in all these works, 
the focus was to study the quantization effects on the CSIT 
under noiseless feedbaclj^. In practice, the CSIT feedback 
may not be error-free due to the feedback channel noise. 
Unlike the forward channel where heavy error correction 
coding can be applied to the time-interleaved payload, the 
limited CSIT feedback has to be received at the transmitter 
with minimum latency and hence, time interleaving is not 
possible. Furthermore, in most systems, the number of bits 
available for feedback is very limited (such as 6 bits in 
WiMAX) and hence, it will be more effective to utilize all the 
limited bits to carry the CSIT rather than wasting some bits 
to protect the CSIT feedback. The issue of noisy feedback 
is considered in Q, For example, in lH, d, 

1 8 1, the authors analyzed the effect of noisy feedback on 
the point-to-point MISO system and broadcast channel @. 
However, the authors did not incorporate the noisy feedback 
into the algorithm design. In ifTOl . ifTTl . the authors design 
a channel optimized quantizer for point-to-point MISO link 
to incorporate the noisy feedback link. However, extension 

'The CSIT feedback index is always received correctly at the transmitter 
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to SDMA system is not-trivial. As illustrated in lfT2l . the 
sensitivity of noisy feedback is much higher in SDMA 
systems and it is critical to take into consideration of noisy 
feedback in the robust limited feedback for SDMA systems. 

When we have noisy CSIT feedback, there may be sig- 
nificant performance degradation because an erroneous CSIT 
feedback will make the base station selecting a wrong precoder 
for the user, which not only decreases the received signal to 
noise ratio (SNR) of the user but also increases the interference 
from other scheduled users since their assigned precoders are 
no longer orthogonal to the target user As we shall illustrate 
in the paper, adopting a naive design approach (design the 
limited feedback codebook assuming error-free feedback and 
testing its performance in the noisy limited feedback situation) 
will result in very poor SDMA performance. In order to obtain 
a feedback-error-resilient design, there are several first-order 
technical challenges to be addressed: 

• Robust SDMA preceding & user selection: Due to 

the noisy limited feedback, the CSIT index received at 
the base station may not be the same as that sent by 
the mobiles. As such, the selected precoder may not 
match the actual CSI at all, resulting in additional spatial 
interference among the selected SDMA users. 

• Robust rate adaptation: To achieve a high system 
goodput advantage, robust rate adaptation is needed to 
control packet errors due to channel outage. 

• Robust Index Assignment: Index assignment refers to 
the mapping of the CSIT feedback indices with the 
precoder entries in the codebook. With noisy feedback, 
CSIT index assignment plays an important role on the 
robustness performance of the SDMA systems. 

• Performance Analysis: Beside robust limited feedback 
designs, it is important to have closed-form performance 
results to obtain useful design insights such as the sensi- 
tivity of CSIT errors in SDMA systems. 

In this paper, we propose a robust noisy limited feedback 
design for SDMA systems. The solution consists of a real-time 
robust SDMA precoding, user selection and rate adaptation 



as well as an offline feedback index assignment algorithm. 
We formulate the robust index assignment problem into a 
Traveling Salesman Problem (TSP) and derive a low complex 
but asymptotically optimal solution. Simulation results show 
that the proposed framework has significant goodput gain 
compared to the traditional naive uncoded design (SDMA 
design assuming error-free and uncoded feedback) and naive 
coded design (SDMA design assuming error-free feedback 
but the limited CSIT feedback bits are protected by FEC) 
under noisy limited feedback channel. Furthermore, we show 
that despite the noisy feedback channels, the average system 
goodput of the proposed robust SDMA design scales as 
^( "nr^-T/'' (C/{) — log2(A^n))) for interference limited scenario 
(Cfb < (jT-T— 1) log2 P+log2 Nn) and 0{nT-log2 P) for noise 
limited scenario where tlt is the transmit antenna number, e 
is the target outage probability, C fb is the number of feedback 
bits, P is the transmission power and Nn is a constant. We 
find that in interference limited scenario, the average system 
goodput increase linearly with the number of feedback bits 
for fixed feedback SER and converge to a constant number 
for fixed feedback SNR. In the noise limited scenario, the 
average system goodput increase linearly with the number of 
transmit antenna and the forward channel SNR (log2 P). 

II. System Model 

In this paper, we shall adopt the following convention. X 
denotes a matrix and x denotes a vector denotes matrix 
hermitian. 

A. Forward MIMO Fading Channel Model 

In this paper, we consider a multi-user system with a base 
station having tit transmit antennas simultaneously transmit- 
ting to riT one antenna active users from a total of K users. We 
shall focus on the case when K > tit so that user scheduling 
in addition to precoder adaptation is important. The base sta- 
tion separates ut data streams to the active users by precoding. 
Each active user k is assigned a n^- x 1 precoding vector 
Wfc. The precoder {wfe}^j!lj^ are a set of unitary orthogonal 
vectors selected from a codebook of multiple sets of unitary 
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influence on it especially in high SNR scenaricQ. Hence in 
this paper, we shall focus on utilizing all the Cfh feedback 
bits on the feedback of channel shape and do not feed 



orthogonal vectors. Let Xk be the transmitted symbol of user k 
with = 1 and yk denote the received symbol of user 

k. The forward channel is modeled as: 



Vk = \ — h 

V riT 



Zk 



(1) 



where P is the transmission poweiQ is the x 1 complex 
channel state vector of the fc*'* user, A is the active user 
set and Zk is the additive white Gaussian noise with zero- 
mean and unit variance. We assume that the transmit antennas 
and users are sufficient separated so that the channel fading 
between different users and different antenna are modeled 
as independent and identically distributed (i.i.d.) complex 
Gaussian process with zero-mean and unit variancj^ 

We consider slow fading channels where the fading is quasi- 
static within a scheduling time slot for each user. This is 
a reaUstic assumption for pedestrian mobility ( 5km/hr) as 
the packet duration is of the order of 500ns (such as Wi-Fi 
and B3G). Due to the quasi static fading and noisy limited 
CSI feedback, there is uncertainty on a user's instantaneous 
mutual information (a function of the instantaneous CSI of 
all K users) at the transmitter. Hence, there exists potential 
packet errors (despite the use of powerful channel coding) 
due to channel outage when the transmitted data rate of user 
k exceeds its instantaneous mutual information. 

B. Limited Feedback Processing at the Mobiles 

In this paper, we consider FDD system and assume the CSI 
is estimated at each user (CSIR) perfectly and fed back to 
the base station through a feedback channel with a limited 
feedback capacity constraint C/b bits per fading block per user. 
The CSIR of user k, h^, consists of two parts: channel gain 



l|hfe|| = yJ2i=iK channel shape hfc = hA;/||hfe||. As 
will be explained in section Hill the average system goodput 
is dominated by channel shape and channel gain has little 

^In this paper we assume equal power allocation since power allocation will 
only bring marginal performance gain under high SNR. This is also assumed 
in (2), (B], [13 etc. 

^In this paper we assume the large scale fading pai'ameters (path loss and 
shadowing) between a base station and all the users in the cell are the same. 



back channel gain ||hfc||. 

We assume the K mobile stations quantize the channel 
shape of the local CSIR with a codebook consists of multiple 
sets of orthonormal vectors: 



M 



(m) 



(2) 



where denotes the quantization codebook with a cardinality 
— N, V'™-' is the m*'* orthonormal set in the codebook 
and M is the number of orthonormal sets which is given by 
M = —. The AI orthonormal bases of are generated ran- 
domly and independently similar to ifTSI^I . Define the distortion 
function between two x 1 vectors Vi , V2 as: 



d(vi,V2) = 1 - |v{v2r = sin^(Z(vi,V2)) 



(3) 



At the k — th mobile, the quantized channel shape hfe can 
be expressed as: 



argmind(v,hfe) 



(4) 



In order to reduce the total feedback overhead for all K 
users, user k shall decide whether to feed back its CSIT or 
not based on the criteria: 



d(hfc,hfc) < S; |hfc||^ > gth 



(5) 



where 6 is the threshold for the distortion from actual channel 
shape to the quantized channel shape and gth is the threshold 
for the channel gain. The motivation of ^ is allocating the 
feedback bits to the users with smaller CSI quantization error 
so as to reduce the potential spatial interference among the 

''when there is only 1 active user in a cell, the system will degrade to 
a MISO system and this claim will be invalid. However, in MISO systems, 
interference is no longer a severe problem and there exists a lot of optimizing 
schemes. In this paper we will focus on strict SDMA system where there are 
more than 1 user served by the BS. 

^Grassmaimian codebook is not a good choice here. Maximizing the 
minimum distortion among M bases is not equivalent to maximizing the 
minimum distortion among all the vectors in the codebook. 
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SDMA Streams. (5 is a system parameter and can be selected 
offline. When the number of total user K is large, 5 can be 
selected as a very small number to decrease the quantization 
error. 

C. CSIT Feedback Channel Model 

Unlike most of the previous literature where the limited 
feedback channel is assumed to be noiseless, we are interested 
in the more realistic case where the feedback channel may be 
noisy. Note that since the CSI feedback has to be delivered in a 
timely manner, effective FEC coding over many CSI feedbacks 
is not possible and hence, the feedback error cannot be ignored 
in practice. 

The CSIT of user k is quantized locally and encoded into 
Cfb bits. The set of CSI indices sent at each of the K mobile 
stations {I^*^'^''} and the corresponding CSI indices received 
at the base station {I^,^'^''} both have cardinaUty of = 
Assume the modulation symbol in the CSIT feedback channel 
has constellation M and the corresponding mapping from 

{ MS) 

the CSIT index 2, to the constellation point m e is 



given by the 1-1 index mapping function A4 — ^(2^^ 



(MS). 



The 



probabilistic relationship between the CSIT feedback symbol 
sent (^(I^*^'^-') and the CSIT feedback symbol received by 
the transmitter ^(I^^"^^) can be characterized by the feedback 



channel transition matrix P 



ch 



{P±,m,}, Where: 



Suppose the condition in (|5]) is satisfied for the k-th mobile, 
the index of hj., t'^^^\ is then mapped to a constellation 
point TOfc using an index mapping function mk — 
and feed back to the base station. Due to the noisy feedback 
channel, selection of the index mapping function ^ becomes 
important and will affect the robustness of the SDMA system. 



III. Base Station Processing: SDMA Precoding, 
User Scheduling and Rate Adaptation 

In this section, we shaU discuss the base station processing 
based on the limited feedback sent from the K mobiles over 
a noisy feedback channel with index transition probability 
PcsiT- Specifically, we shall discuss the SDMA precoding, 
user selection and rate adaptation. 

A. System Goodput 

Consider a full multiplexing system where base station 
schedules ut active users for SDMA. Define A as the active 
user set and as the precoder for a scheduled user k. We can 
write the instantaneous goodput (b/s/Hz successfully received 
by user k) pk for user k as: 

Pk - Tk (if ) • 1 [rfc (if ) < Cfc (hfc , Wfc )] (8) 
where is the data rate of the packet of user k and is a 



pch ^ p 



e(lf^))-mne(lf^^)=™fc 



function of received CSIT at base station, Ck is user ^'s 

,mi,mk e M. ,(6) . ^ , , . 

instantaneous mutual information and 1(A) is an indicator 



Hence, the reliability of the CSIT feedback channel is char- 
acterized by the CSIT feedback channel transition matrix 
Pch- Note that Pch depends on the average feedback SNR, 
feedback constellation and so on and can be offline evaluated 
analytically or numerically through simulations. 

Given Pch and the index mapping rule ^, the stochastic 
relationship between the CSIT index sent by the mobile station 



function which is 1 if the event A is true and otherwise. 
User k's instantaneous mutual information Ck is a function 
of its channel state as well as the assigned precoder w^. 
Specifically, Ck can be written as: 



Cfc(hfc, Wfc) = log2(l 



) (9) 



t 



(MS) 
k 

r{BS) 



and the CSIT index received by the base station 
is characterized by the CSIT index transition probability 



Pcs/T = K%^J given by: 



pCSIT ^ 



riBS) _ 



(MS) ^ ■ 



(7) 



We define = Z(hfc,Wfc) as the angle between the actual 
channel shape and the assigned precoder for user k. Similarly, 
define = Z(hfc, V2-(ms)) and ip{vi,\-j) — Z{vi,Vj) as the 
angle between the actual channel shape and the quantized 
vector of user k and the angle between two x 1 vectors 
V,:, V,. 
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Consider high effective SNR asymptotic scenario where 
;^||hfep is sufficiently large. We can further simplify (|9]l to: 



Cfc 



l0g2(l + 



(10) 



Sj5:tfc,je>tCOs2(p(hfe,Wj) 
As shown in Figure [T] the approximation is quite good for 
moderate to high SNR. 

Remark 1: While the approximation in ( fTOb fails when 
J2j^k jeA '^^^^ '^'(hfc, Wj) — > 0, it will not affect our design 
and analysis because for a practical target PER (e.g. 10^^)^ 
when J2j^k jeA'^'^^^ ^ 0, both Ck in ^ and 

([Tol l will be large enough and no outage will occur In other 
words, the case J2j^k jeA '^'-'^^ "y^C^fe, Wj) — > will not be the 
performance bottleneck and there is no loss of generality to 
focus on the bottleneck case when J2j^k j£A '^'-'^^ fO^kj'vvj) 
is not close to 0. 

Remark 2: From equation (fTol l. we can see that at high SNR, 
the instantaneous mutual information does not depend on the 
channel gain ||hfc|| and hence the average system goodput at 
high SNR is dominated by channel shape. 

The average goodput for a scheduled user k can be ex- 
pressed as: 

Pk = Kj-iBS)Eh^[pk] 

= E^(«.)[rfe(lf^)).Pr(rfc <Cfc|lf^))] (11) 

In practice, there is a target PER requirement e associated 
with different application streams (e.g. e = 0.01 for voice 
applications), which is expressed as: 



l-Pr{rk{li'''^)<Ck\ir')^^ 



(12) 



The average system goodput for all scheduled users can be 
written as: 

keA 

B. SDMA Precoding and User Selection 

Generally speaking, based on the CSIT feedback, the active 
user set A and corresponding precoder shall be jointly 
optimized to maximize the overall average system goodput in 
(fT3T l. Assume that the number of total user K is large enough 



such that the base station can always fuUy schedule active 
users. A natural algorithm includes an exhaustive search over 
all the combinations of riT users out of a pool of K users 
and jointly optimize the riT precoders for the selected users 
to maximize the average system goodput. However, exhaustive 
search has exponential complexity in terms of number of users 
and is not practical as an online algorithm especially for large 
number of users. Therefore we shall adopt a simple orthogonal 
scheduling and precoding algorithm. At base station, take the 
received CSIT v (bs) as the estimation of h^, the term cos^ 6 
in the signal power term in (fTOl i is maximized if — v ibs) . 
Similarly, the interference term J2j^k jeA'^'-'^'^ f0^k,'v^j) is 
minimized by choosing w, _L v (ss) , j ^ k,j e A. 

It can be easily shown that the above two equations can be 
satisfied simultaneously with: 



v„(BS) -L V, 



.(BS), 



Vj ^ k,j,ke A. 



(14) 



In other words, the received CSIT of tit scheduled users 
{v (BS)} form an orthonormal set V*^™-* in the codebook J^. 
In fact, SDMA with orthogonal precoding is also called per 
unitary basis stream user and rate control (PU2RC) II16I and 
has been widely used in the standards such as 3GPP-LTE[3]. 
The main feature of PU2RC is that it could accommodate 
limited CSIT feedback in a natural way. For instance, the 
multiuser precoders are selected from a codebook of multiple 
orthonormal bases. The importance of PU2RC for the next- 
generation wireless communication motivates the investigation 
of its performance in this paper. In this paper, we consider a 
simplified PU2RC system where scheduled users have single 
data streams, which are separated by orthogonal precoders. 

The precoding and user scheduling strategy is briefly sum- 
marized as follows: 

• Step 1: search N/tit groups of orthonormal sets V*-™^ in 
for one set V'™-** in which each vector is received by 
the base station as a CSIT feedback for at least one user 
If there exists multiple satisfying sets, randomly pick one. 

m Step 2: For an orthonormal vector Wk G V*- '"-'*, randomly 
select one user from the group of users with received CSIT 
Vfc as the scheduled user and set as its precoder 
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• Step 3: Repeat Step 2 for all the vectors in V^™-'* 

C. Robust Rate Adaptation 

With the proposed preceding and user scheduHng strategy, 
Cfe in ( [Tol l is given by the following lemma: 

Lemma 1 (Mutual Information at High SNR). At high 
downlink SNR, when the scheduled users are using nx 
orthonormal precoders to transmit, the mutual information of 
a scheduled user is given by: 



Ck = -21og2(sin( 



Note that in practice, i* and Ns'^ {^i^^^ ) can be offline pre- 
calculated from the channel transition matrix p'^^^'^ and the 
distortion between the codewords sin ipi,j. For example, sup- 
pose our codebook is given by F = {vi, V2, V3, V4}. Assume 
e = 0.1 and Pg^^^ = 0.70, sinZ(vi, vi) = 0;P^i'^^^ = 
0.10,sinZ(vi,V2) = 0.5; Pg^''^ = 0.11, sinZ(vi, V3) = 
OA;Pg^^^ = 0.09, sin Z(vi,V4) = 1, then Ns"-^{1) = 
{v2, V3} and i* ~ 3. 

Using Lemma ID and define e^es = ( — 

E-^TTT P'^^/^q\, the transmission rate of user k is 



(15) given by: 



Proof 1. Please refer to Appendix A for details. 

Tk can be calculated from the requirement of the conditional 
PER target in (fT2] l. which is given by: 

Pout = Pr(-21og2(sin0) < rfc|lf = e (16) 

Yet, one critical challenge in solving for ru in ( fTSI l is the 
knowledge of the CDF of siwO conditioned on the CSIT 
■j^iBS) ^ This is in contrast with the conventional approach 
of maximizing the ergodic capacity in which only the first 
order moment of the random variable sin 9 is needed. In the 
following lemma, we shall give a tight upper bound on the 
conditional CDF of sin 6*, which is critical to solving for a 
closed-form rate adaptation solution. 



rk 



pCSIT 



)2(,.r-i) ^ sin ip (BS) 



(18) 



IV. Robust Index Assignment 

Index mapping algorithm is important when there is noise 
on feedback channel. In this section, we shall optimize the 
index assignment function ^ to maximize the system goodput 
in ( fT3] ). This is equivalent to minimize the feedback distortion 
between the channel shape hk and corresponding precoder 
of user k which is sin^ 6. This feedback distortion is 

contributed by two parts: distortion from quantization sin^{(j)) 
and distortion from feedback error sin^((^(v^^'^'', v^*^"^-*)). 
Based on: 



Lemma 2 (Upper bound of the PER Pout)- The conditional 
PER Pout of the forward channel ( conditioned on the limited 
CSIT feedback received at the BS) is given by: 



(2-* -sin(p 

TjCSIT 



CSIT 



BS) 



j&Ns' IBS) 



(17) 



where i* = argmax^^j^^^^^^BS).^ sinip^(BS) ^ and iVs'^(I^ ') 
is the set of neighboring codewords of l^^^^ satisfying: 



Proof 2. Please refer to Appendix B for details. 



(Pj-(BS) j.{MS) — (j) < < Lp^iBS) ^{MS) + 



(19) 



and sin cj) < S where 6 is chosen to be a small value to 
avoid excessive spatial interference among the SDMA streams. 
As a result, we shall omit quantization distortion sin^ (j) and 
focus on minimizing the distortion introduced by feedback 
enov d{-v[^^\w[^'^^'') = sin^{ip{v\f^\w[^'^^'')). The average 
distortion is given by: 



N N 



t=i j=i 



d(v„vj) (20) 



Searching for optimal index mapping function ^ can be 
summarized into the following problem: 
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Problem 1 (Robust Index Assignment Problem). Find an 
optimal index assignment function to minimize the average 
distortion introduced by feedback error: 

N N 

C - argmin^^Pr(vO • P^^^^^^^-^ ■ d(v„v,) (21) 

1=1 j=i 

where Pr{-Vi) is the probability that is the quantization 
output for hk and d{'Vi,'Vj) is the distortion between two 
codeword given in 

Note that the insight of the above formulation is that a good 
index mapping function shall map 2 precoders vi,V2 G J- 
with smaller distortion d{vi,V2) to the constellation points 
mi , 1712 with larger transition probability P.^!^_^ .^^^ . 

In general, finding the optimal mapping ^(.) involves 
combinatorial search. When the number of feedback bits is 
small, the computation complexity of exhaustive search is still 
acceptable. However, when the number of transmit and receive 
antennas gets larger and more feedback bits are required, the 
exhaustive searching time will increase double exponentially 
with Cfb- This motivates the study on the low-complexity 
solution of the problem. 

Consider a special case when the CSIT feedback index 
is modulated by one iV-PSK symbol. When feedback error 
occurs, the erroneous symbol is likely to be one of the adjacent 
neighbors of the feedback A^-PSK symbol, which is referred 
as the nearest constellation error. The average distortion intro- 
duced by feedback error in (ISTT i can be simplified to: 

,^^(ve-i(i))'^(ve-i(i)'V«-i(2)) , 



m) = Pe-{i- 



Pr(v^-i(2))d(v^-i(2),v^-i(i)) 



) + ••• + 



Fr(v^-l(A,))d(v^-l(Ar),V5-l(l)) 



)} 



(22) 



where P^ denotes the symbol error rate (SER) of the feedback 
channel and D{(_) = K{d) is the average distortion. Assume 
that {Ndi, ■ ■ ■ , Nd^} are N virtual cities, and the distance 
between the virtual cities Ndi and Ndj is given by: 

Pr(vi)d(vj, Vj) + Pr(vj)d(vj, Vj) 



Dis{Nd,,Ndj) = 



(23) 



Equation ( |22] ) can be expressed in terms of distance between 
virtual cities as follows. 

D{0 = i?^s(iVd^-l(Jv),A^dc-l(l)) + 

Dis{Nd^-i(i),Nd^-^2)) + ■■■ + 
Dis{Nd^-i(^N-i),Nd^-i^N)) (24) 

Hence, the optimization metric is equivalent to the total 
distance of a Hamiltonian cycle |17|. From equation ( l24l l. the 
index mapping problem in Problem[T]is equivalent to searching 
shortest path in a Hamiltonian cycle and this can be cast into a 
traveling salesman problem (TSP). This is summarized below. 

Problem 2 (Traveling Salesman Problem). Given a number 
of cities {Ndi, Nd2, ■ ■ • ,Nd]\j}, and the costs of traveling 
from any city to any other city {Dis{Ndi, Ndj)}, what is 
the round-trip route [^^^(1), ^^^(2), • • • ,C~^(-/V)] that visits 
each city exactly once and then returns to the starting city to 
minimize the total distance i\24\l . 

TSP is found to be an NP-hard (nondeterministic polyno- 
mial time) problem and yet, there are a number of efficient 
searching algorithms for the TSP such as the cutting-plane 
method [TSl and genetic algorithm lfT9l . 

In this paper, we propose a simple construction algorithm, 
namely the Circled Nearest Neighbor Algorithm(CNNA). The 
CNNA algorithm is described below: 

• Step 1: Start the TSP travel from a randomly selected 
node Ndi. 

• Step 2: Go to the nearest unvisited node from Ndi. 
If there exists more than one such nodes, we select 
the one with smallest sum distortion to the previously 
visited nodes. For example, if S is the set of nodes 
already visited, and Afi is the set of unvisited neighboring 
nodes of Ndi, we shall select Ndj via the criteria: 
Ndj = argmin^vdjeM J^Ndkes dis{Ndk,Ndj) 

• Step 3: Repeat Step 2 till all the nodes are visited. Then 
go back to the start node. 

Suppose we require that all feedback error results in the 
"nearest neighboring precoder" and we would like to find an 
index assignment such that this can be realized for "2-nearest 
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neighbor error feedback channel". To do that, we assume all 
the neighboring regions of codewords are "equi-probable" and 
isotropically distributed on the surface of a unitary hypersphere 
(hfe-space). For the given topology of the partition region, the 
above algorithm of index assignment will result in a circled 
pattern as illustrated in figure |2] and is similar to the following 
analogy. Suppose we start from north pole of the Earth and 
travel around the world along the latitudes to the south pole. 
When we finish traveling along a latitude, we go down to the 
next one until we reach the south pole. After south pole is 
arrived, we go back to the north pole directly. Define dmin as 
the minimum distortion between two precoder, we have the 
following lemma: 

Lemma 3 (Asymptotic Optimahty of CNNA Algorithm). 
For N-PSK constellation with nearest-constellation error ap- 
proximation, the index mapping solution given by the CNNA 
algorithm ^* is asymptotically optimal for sufficiently large N. 
i.e. 

lim T-;— ^ drain- (25) 

Proof 3. We provide a sketch of proof due to page limit. 
With the proposed algorithm, all error events shall result in 
nearest neighbor codeword errors (for 2-nearest constellation 
feedback channel) except for the codewords serving as the 
starting and ending nodes. Hence we have 

D = N ■ drmn + C (26) 

where c is a constant. When N is sufficiently large, the average 
distortion of the travel is limAr^+oo ^ ~ dmin. Hence the 
proposed algorithm is asymptotically optimal. 

V. Performance Analysis 

In this section, we shall focus on obtaining the asymptotic 
goodput performance under noisy limited feedback. The trans- 
mission rate for a scheduled user in ( fTSl l can be bounded by: 

-21og2(55 + sin(/7^(f5S) ..) <rk< ~2\og2{s\nip^{BS) .,). 



In fact, the upper bound and lower bound of above 
are both very tight when 5 is small. Asymptotically, the 
two bounds will meet each other as the number of feedback 
bits Cfb goes to infinity since 6 will approach to and 
sin(^_(i3S) is the dominant factor to rfc. Since we are 
interested in the first-order analysis, the data rate Vk can be 
taken as: 

Tfc = 0(-21og2(sin(^ (BS) )) (27) 

k ' 

Substitute ( |27] | into (fTsl i. the average system goodput is 
given by: 

g = 0(-2nT(l-e)log2(E (Bs,[sin(p (i,s),J))(28) 

From ( l28T l. we can see that system goodput performance 
depends on the worst case distortion sin 1^9 (bs) of an index 
assignment function ^. Define N^ii) as the set of neighbor- 
ing points of a constellation point i (including itself) with 
^j£N'{i)^<^h{iTi) = 1 — e. The cardinality of the set is 
Nn = \Nr,{i)\- Note that N^ii) is the Nn largest terms in 
the ith row of Pc/iS As a result, both P^. and Nn are two 
first order parameters to characterize the quality of feedback 
channel. We have the following lemma: 

Lemma 4 (Lower bound for E_(bs) (sin(p.j-(BS) ,)). For suf- 

■^k -^k '* 

ficiently large N, we have: 

E (BS)(sinv^ (BS) > (^)^^^ (29) 

k k ' iV 

Proof 4. Please refer to Appendix C for details. 

Numerical results show that with optimal or near optimal 
index assignment, the lower bound is quite tight as illustrated 
in figure [3] Take the lower bound in ( |29] l as an approximation 
of E (BS) (sin (/3 (BS) and substitute into ( |28] |. we have the 

-^k ^k 

following theorem: 

Theorem 1 (Asymptotic System Goodput Q in Interference 
Dominant Scenario). When the number of feedback bits Cfb 
satisfies Cfb < {tit — 1) log2 -P + log2 -^n, the system is 

"Example: For 8PSK with lOdB feedback SNR we have Nn = 3 when 
e = 0.03. 
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dominated by interference. For sufficiently large transmission 
power P and quantization codebook size N, the average system 
goodput is given by: 

g = o(!^lii_!)(C/,„log2(7V„))) (30) 

Corollory 1. In interference dominant system with sufficiently 
large transmission power P and quantization codebook size N: 

• With fixed SER Pg on the feedback channel, Nn is a finite 
constant and the average system goodput is given by: 

g = oC^i^c,,) (31) 

• With fixed feedback SNR, Nn scales with N as iV„ ~ 
c ■ where c is a constant and the average system 
goodput is given by: 

Q = OC^I^) (32) 

riT — 1 

The goodput order of growth results in ( [3T1 ) and ( [32] ) are 
also verified against simulations in Figure |4] and Figure |5] 

On the other hand, when the number of feedback bits is 
sufficiently large, the SDMA system will operate in the noise- 
dominated regime. 

Theorem 2 (Asymptotic System Goodput Q in Noise Domi- 
nant Scenario). When the number of feedback bits C fb satisfies 
Cfb > {riT — 1) log2 P + log2 Nn, the system is noise- 
dominated. For sufficiently large transmission power P and 
quantization codebook size N, the system goodput is given by: 

g = Oinr-log^P) (33) 

Proof 5. Please refer to Appendix D for details. 

The above result also reduces to that for SDMA system with 
noiseless feedback when iV„ = . Note that Nn ^ 0{1) 
actually corresponds to an asymptotically noiseless feedback 
channel. We shall note that Cfb > {nr ~ 1) log2 P + log2 Nn 
may not be able to be satisfied by simply increasing Cfb 
since Cfb will be canceled when Nn = c ■ 2^f'' (e.g. constant 
feedback SNR). This indicates that one could not enhance the 



CSIT quality by increasing Cfb if the feedback SNR is kept 
constant. 

VI. Results and Discussions 

In this section, we study the performance of the proposed 
robust SDMA system under noisy limited feedback. We com- 
pare the performance of proposed design with two naive 
designs under the same feedback cost. In the uncoded naive 
design, the SDMA system is designed as if the limited CSIT 
feedback were noiseless and the limited CSIT feedback bits 
are uncoded. In the coded naive design, the SDMA system 
is similar to the uncoded naive design except that the limited 
CSIT feedback bits are protected by hamming code. In the 
simulations, we consider an SDMA system with = 4 
forward SNR 20dB and K = 100. We set the thresholds 
(5 = 0.1 and gth = 2. 

A. System Performance with respect to the Number of Feed- 
back Bits Cfb 

Figure |4] illustrates the average goodput versus the number 
of feedback bits of the SDMA system with a fixed feedback 
SER of 0.^ There are significant goodput gain compared 
with both the naive uncoded and coded designs. The order- 
of-growth expressions in Theorem [T] and Corollary [T] are also 
verified. 

Figure |5] illustrates the average goodput versus the modula- 
tion level per feedback symbol with a fixed average feedback 
SNR of lOdB and a fixed number of feedback symbols. There 
are also significant goodput gain in the proposed scheme 
compared with both the naive uncoded and coded designs. 
In addition, there is a tradeoff relationship in the feedback 
constellation level for the baseline systems. With lower feed- 
back constellation level (such as BPSK), the feedback is more 
robust but the average goodput performance is limited by the 
resolution in the CSI feedback. On the other hand, for large 
feedback constellation, the average goodput performance of 
the reference baselines are poor because the performance is 
limited by the feedback error 

^This corresponds to a feedback SNR of 10 dB when Cfi, = 3 and 8PSK 
is adopted. 
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B. System Performance with respect to Feedback Quality 

The feedback quality can be specified by feedback SER and 
feedback SNR. Figure |6] shows the average system goodput 
versus SNR with different feedback SER and fixed feedback 
bits Cfb — 8. It is shown that with the proposed design, the 
system goodput decrease much slower with the increasing of 
the feedback SER. In figure |71 we show the average system 
goodput versus feedback SNR with fixed feedback bits Cfi, — 
6. With proposed design the system goodput increases much 
faster with the increasing of the feedback SNR. 

VII. Conclusion 

In this paper, we proposed a robust noisy limited feedback 
design with a joint user scheduling and precoder scheme as 
well as rate adaptation and robust index assignment optimiza- 
tion algorithms. We convert the index assignment optimization 
problem to a Traveling Salesman Problem(TSP). Simulation 
results show that the proposed framework has significant 
goodput gain compared to the uncoded and coded naive 
designs. Furthermore, we show that despite the noisy feed- 
back, the average system goodput scales as O ( "^^"'^j^ (C/;, — 
log2(iV„))) and 0{nT ■ \0g2P) in the interference limited 
regime {Cfb < {ut — 1) log2 P + log2 Nn) and noise-limited 
regime respectively. 

Appendix-A: Proof of Lemma[T] 
From equation (|9|l, the instantaneous mutual information can 



be written as: 

Ck{h.k,wk) w log2(l 



cos^e* 



-) 



= -21og2(sin0) (34) 

(l34l i is because all the scheduled precoders Wj forms 
an orthogonal bases of the ut dimensional space and 
(cos (p{h.k,Wji), ■ ■ ■ , cos iy9(hfc, Wj„y )) is a vector on the unit 
sphere of the ut dimensional space, where w^i to Wj„y are 
fiT scheduled orthonomal precoders. Hence we have: 



008^6'+ cosV(hfc,Wj)) = 1. (35) 



Appendix-B: Proof of Lemma|2] 

Given received CSIT Ij. , the outage probability can be 
written as: 

Pr(2^"2" — simp (BS) < sin cj) < ^/6) 



Pr(sin (j) < VS) 



■P 



i&Ns- IBS) 



CSIT 



(36) 



where i* = argmax.^^^^^^jf-ss)^ siiiip^(BS) . and Ns'^{ll,^^^) 

' riBS) 



is the set of neighboring codewords of X^. satisfying: 

y~^-^7\r ^ /-T-iBS)^ P^'^/jK-^ ^ 1 — 6 and ^ f^(BS)^ P^^/f^\ — 

P^/\ls) < 1 - £■ From EOj, we have the CDF for sin cj): 



Pr(sin0 < a;) = a;2(«r-i) 



(37) 



Subscribing ( |37] | into (l36T l, we can simplify the outage prob- 
ability as: 



Pout < (1 



(2-* -sin(p.^(BS),j2(nT-i) 



E 



J(nr-l) 
TjCSIT 



jeNs' IBS) 
k 



Appendix-C: Proof of Lemma|4] 

The asymptotic expression for E (bs) (sin (p^(BS) ) is de- 
duced assuming all the channel shape quantized to the ele- 
ments in A^s'^(Ij, ) forms a neighboring region of v^(b 



S) 



defined as a hypersphere: 



Cbiv (BS)) {hfc : sinZ(hfc,v (BS)) < rcj (38) 

k k 

where rc^ is the radius of the hypersphere and v^. lies at 
the edge of the hypersphere. Since Pr(sin Z(hfe, < x) — 
^2(nT-i)^ and the size of the hypersphere is: 

Pr(sinZ(h.,v^,B.))<rcJ = ^^^^^^|^>^ (39) 

where |iVs^(I^^^^)| is the cardinality of Ns'{ll^'^'^). There- 
fore we have rc^ > 



2(,^J,-l) 



Since v^. lies at the edge of the hypersphere, we have: 

E (BS)(sin(/? (BS) ) > (-rf )'<"r-'' . 

-^k -^k iV 
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Appendix-D: Proof of Theorem |2] 

Consider the situation with noise power larger than in- 
terference power 1 > E[;^||hfe|p sin^ as noise dominant 
scenario. With optimal or near optimal index assignment, 
we have Elsin^ 9] < (^)^' 



^ — Instantaneous Mutual Information 
+ Approximation 



Therefore, we require P 



-j^j . -1 < 1 which gives C/b > (riT-l) logj P + iog^ Nn. 
The instantaneous mutual information of a user in (|9]l can 



-I 



hkP COS^ I 



As the 



be simplified into Ck ~ log2(l + 
channel gain ||hfc|p scale with nx, without loss of generality, 
we can also select gth in dU at order 0{nT)- Therefore the 
transmission rate for user k has the order ^ 0(log2 nr) 
hence the system goodput Q has the order: 



g ~ 0(nT -logaP). (40) 
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Goodput vs Feedback Bits — Fixed Feedback SER 




An Example oF 
Index Mapping 



Simulated Goodput with proposed design (NPSK)-^ 
Simulated Goodput with proposed design {NQAM} 
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Order-of-growth expressions . 
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Naive Coded Design 



Naive Uncoded Design 
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Number of Feedbacl< bits C,^ 



Fig. 4. System goodput versus number of feedback bits. Naive Uncoded 
Design refers to design assuming the limited feedback channel is noiseless 
and the feedback bits are uncoded and Naive Coded Design is similar to naive 
uncoded design except the feedback bits are protected by Hamming code. The 
SER on feedback channel is fixed at 0.2. The dotted curve of naive coded 
design is not smooth because of the overhead in the parity bits of Hamming 
code {k, 2*) used to protect the CSIT feedback. 



Goodput vs Feedback Bits — Fixed Feedback SNR 



Fig. 2. A heuristic TSP travel on a hypersphere. 
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Fig. 3. Accuracy of approximation on /(^) 
Lemma [4] 



E^-Cfss) {sinip IBS) in 
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Feedback Consteliation 



Fig. 5. Average system goodput versus Limited Feedback constellation 
Level. Naive Uncoded Design refers to the limited feedback design (assuming 
noiseless feedback) and the feedback bits are uncoded and Naive Coded 
Design is similar to naive uncoded design except the feedback bits are 
protected by Hamming code. In all the schemes, the feedback SNR(lOdB) 
and the number of feedback symbols are fixed. The dotted curve of naive 
coded design drops dramatically when the number of feedback bits gets large 
because Hamming code only protects nearest neighbor errors. 
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Goodput vs SNR (SER = 0.01 ,0.2. 0.5 respectively) 



- Proposed Design, SER = 0.01 
Naive Uncoded Design, SER = 0.01 
Proposed Design, SER = 0.2 
Naive Uncoded Design, SER = 0.2 
Proposed Design, SER = 0.5 
Naive Uncoded Design, SER - 0.5 
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Fig. 6. System goodput versus Forward SNR at different level of feedback 
SER. Number of feedback bits Cfi, = 8. With proposed design, the system 
goodput decreases much slower with the increasing of feedback SER. 
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Fig. 7. System goodput versus feedback SNR. Number of feedback bits 
Cfi, = G and forward transmission SNR is 20 dB. The performance of 
the naive design is very bad since the Nn precoders that may be adopted 
when feedback en'or occur ai'e randomly selected and the transmission rate is 
determined by the worst precoder. 
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